On the emergence of critical behavior in evolving financial networks 
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We introduce a simple model for addressing the controversy in the study of financial systems, sometimes taken 
as brownian-like processes and other as critical systems with fluctuations of arbitrary magnitude. The model 
considers a collection of economical agents which establish trade connections among them according to basic 
economical principles properly translated into physical properties and interaction. With our model we are able 
to reproduce the evolution of macroscopic quantities (indices) and to correctly retrieve the common exponent 
value characterizing several indices in financial markets, relating it to the underlying topology of connections. 
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The introduction of the theory for brownian motion in the 
analysis of financial data in 1973 by Black, Scholes and 
MertonO]] boosted a new field of research in economy, ex- 
plaining the price evolution in an organized market. For that, 
an analogy between the fluctuating forces that a Brownian par- 
ticle experiences in its interaction with the surrounding en- 
vironment was introduced corroborating previous studies on 
portfolio management^ and establishing a theoretical frame- 
work that was extendend from interest ratelst] to credit riskG ] 
evaluation. More recently, a procedure was introduced 10, li] 
for quantitatively describing financial data, namely by explic- 
itly deriving a Fokker-Planck equation for the empirical prob- 
ability distributions which catches the typical non-Gaussian 
heavy tails of financial time-series. Such description sup- 
ported a cascade model for the probability density function 
deformation accross scales. However, as Mandelbrodt al- 
ready pointed out in the 60's[7j|, other typical features are 
observed in the variation of prices, namely scale invariance 
behavior, which cannot be explained by means of a cascade 
model. Therefore, the modelling of financial index dynam- 
ics and other complex phenomena taking into account non- 
Gaussianity and scale invariant behavior has been indicated 
as an appealing problem to address in the scope of finance 
analysis and statistical physics JH. 

In this Letter we propose a model based on first principles 
borrowed from Economy which are straightforwardly trans- 
lated to the physical context, dropping the assumption for the 
system to be in equilibrium, differently from the Brownian 
particle approach. Since the main concern in risk manage- 
ment deals with the distribution of drops on financial market 
indices, we restrict ourselfs to that distribution. We show that 
using such first principles, we are able to retrieve the typical 
behaviour of financial indices drops and furthermore to derive 
quantitatively the exponent characterizing the critical behav- 
ior in the system. We also show that this exponent depends 
linearly on the underlying complex network exponent. 

Our model is based in the paradigm of self-organized crit- 
icality (SOC)|@] where each local avalanche represents one 
local economical crash. Such local crashes appear at no par- 
ticular moment and may cascade into an avalanche of arbitrary 



size, i.e. a succession of an arbitrary large number of local 
economical crashes, taken as a global economical crisis. 

The model considers a collection of agents, labelled i and 
j in Fig. Q] which establish among them bi-directed connec- 
tions. When agent i delivers labor Wi to agent j (labor con- 
nection) it receives in return a payment proportional that labor, 
Ej = a.ijWi (energy connection). Since the payment Ej is in 
itself a part of the labor produced by agent j, both Ej and Wi 
have units of energy. 

The factor is a (adimensional) measure of the labor 
price, defined as 
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where k ou t,i and ki n j are the outgoing and incoming 'la- 
bor' connections of agent i respectively. With such defini- 
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FIG. 1: Left: Illustration of an economical connection between two 
agents i and j. Agent i transfer labor Wi to agent j receiving an 
energy Ej = oujWi where ouj measures how well the labor is re- 
warded. This interaction attributes to agent i an amount of "internal 
energy" Ui = Wi — Ej that can be accumulated by the agent up to 
a threshold Uth beyond which it is distributed among the neighbors 
(see text). Right: The transfer of labor is done according to a prefer- 
ential attachment scheme: the agent prefers to work for agents which 
have already a large number of labor connections to them. 
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tion otij G [0,2], where the middle value o>ij = 1 occurs 
when the labor is payed with the same amount of energy. For 
> 1 the labor of agent i is payed above the amount of en- 
ergy Wi it delivers, i.e. agent i profits from the trade, whereas 
for ctij < 1 the opposite occurs. According to basic econ- 
omy principles 1IT2I1 such over- and underpayments are ruled 
according to the demand of and supply from a specific agent, 
which is properly incorporated in our definition in Eq. (Q3: a 
large (small) ki n j indicates a large (small) supply for agent i 
and a large (small) k ou t,i indicates a large (small) demand of 
agent i. Thus, the difference k out .i — fcj n j measures the bal- 
ance between the demand of an agent with the supply it has. 
For the two limit cases of a very large supply and a very large 
demand, i.e. k in ,j > k ou t,i and fc mj <C k outyi , the price a tj 
approaches the two limiting values, and 2 respectively. 

At each agent i, the energy balance implies the introduction 
of an additional internal energy defined as U = J2out — 
Ej) + J2in ii^i ~ Wj). Here, we assume a mean-field ap- 
proach where all a^ in the system are substituted by the aver- 
age value a = (otij) yielding for each agent 



Ui = (1 - a){k oll 
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These rules are incorporated in a system of N agents, ini- 
tially with the agents having a fixed number k out of out- 
going connections according to a preferential attachment 
scheme 111 611 . Initially the outgoing connections follow a 6- 
distribution P ou t(k) = S(k — k out ) and the incoming con- 
nections follow a scale-free distribution Pi„(k) = k r . 
Behind this topology underlies the empirical observation in 
economic-like systems that agents are more likely to deliver 
their work to agents receiving already significantly amount of 
work 117 1 . 

When the system evolves, the number of agents remains 
constant but at each event-time n one new connection joining 
two agents is introduced. Both agents are independently cho- 
sen according to preferential attachment schemes as above. 
Consequently, the joint distribution of outgoing and incoming 
connections follows P(k out ,k ln ) = k~J t ° ut k~? in . 

In a single trade, the internal energy of an agent i is ei- 
ther Wi — Ej or E% — Wj depending if the agent is deliv- 
ering or receiving labor. Thus, for m < (resp. Ui > 0) 
the agent increases its asset (resp. liability). Since new con- 
nections are being introduced, the price values (ay) change 
during the system's evolution, eventually promoting the accu- 
mulation of debt in a specific agent. Such accumulation can 
only exist up to a certain threshold Uth,i that depends on how 
much debt th e ag ent can afford. The ability to accumulate 
debt is related ! 1811 to how much influence the agent has in the 
economy. Such influence can be measured by the turnover 
Ti = kout,, + hn,i- Thus, we fixe a threshold d th = Uth,i/T it 
assumed constant for each agent, and during evolution keep 
track of the quantity di = Uj/Ti, we call percentual deficit. 
One sees that the smaller the percentual deficit of one agent 
the larger its consumption in comparison to what it produces. 

Under these assumption we consider that when di < dth the 
agent collapses and an avalanche takes place. This collapse 
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FIG. 2: (a) Evolution of the variation of the total internal energy, 
which shows (b) probability density function and (c) avalanche size 
distribution similar to the ones observed for empirical data (compare 
with Fig.fTJ. 



induces the removal of all the ki n i consumption connections 
of agent i from the system. Consequently, 



U ->■ Ui + (1 - a)h 
Ui 



Uj - (1 - a) 



Ti 



-> Tj - 1, 
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where j labels each neighbor of agent i. If the avalanche 
stops here one has an avalanche of size s = 1. If the collapse 
of the first agent induces neighboring agents to collapse either 
the size of the avalanche increases. 

Figured shows a sketch of the evolution of a typical time- 
series for the logaritmic returns of the total internal energy, 
dU/U. As can be seen from Fig. [2}} the distribution of the log- 
arithmic returns is non-Gaussian with the heavy tails observed 
in empirical dataj^]. In Fig. [2fc the cummulative distribution 
A(s) of the avalanche size s is plotted showing a power-law 
whose fit yields A(s) ~ s~ m with an exponent m = 2.51 
(R 2 = 0.99) which does not change for systems with different 
number of agents. To be comparable to empirical series (see 
below) we consider in our analysis a sampling of data which 
takes one measure of the original series from the system each 
five iterations. 

The condition for the system to remain at the critical state 
can be derived as follows. Using the method of Barabasi- 
Albert-Jeongljjl for the iterative procedure of preferential at- 
tachment, we can derive the degree distribution as P(k) ~ 
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FIG. 3: (a) Time evolution of the logaritmic returns of the DJIA in- 
dex (partial). The box is zoomed in (b) to emphasize the constrast 
between real data (solid line) throught time t and the succession of 
events (dashed line) where the variation changes sign, (c) Probabil- 
ity density functions for some important financial indices, including 
interest rate options (CBOE). (d) Distribution of avalanche sizes de- 
tected throught the evolution of each financial index showing critical 
behavior with (e) an exponent m approximately invariant and simi- 
lar to the one obtained for branching processes and to our model (see 
tex). 



Kok~ 7 , where k is the degree, 7 the exponent of the de- 
gree distribution and Kq is the initial number of outgoing 
connections a node has (here Kq = 1) valid for both the 
consumption and production networks. Since the system is 
in a critical state whenever the expected number of neigh- 
bors that experience collapse (di < d t h) is one or larger, one 
needs to derive the probability for a neighboring node to col- 
lapse. From Eqs. ([3]) a collapsing node i has all its consump- 
tion connections removed from its ki n ,i neighbors, produc- 
ing a reduction in the total number k out ^ of production con- 
nections for its neighbors. Thus, the collapse of one neigh- 
boring node fulfills k out .i - k in .i > d t h(k ut.i + k in ,i) and 
k ut,i - 1 - k in ,i < d t h{kout,i - 1 + k in .i), yielding 

uk in .i < k outil < wki n> i + 1 (4) 

with u> = jj^j ■ Taking a collapsing node, the probabil- 
ity Pbr for a neighbor to also collapse is the probability for 
the above condition to be fulfilled, i.e. Pb r = P{u>ki n ^ < 
k ut,i < wfcin,i + 1), and since k ou t.% is an integer one obtains 
Pbr ~ K^ujkin.i)' 1 ■ Therefore, assuming that the system 
cannot consume an infinite quantity of energy, the expected 



number of collapsing nodes from a starting one is given by 
the condition 



^ kin,iP (kin,i) Pbr — 1 



(5) 



which yields Y^k in i= \ kin.i = (t?o) ' *' e ' ^ e s y stem re ~ 



mains in the critical state lll5ll as 
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where £ is the Riemann zeta-function. Condition ((6) close our 
model, since uj depends only on d t h and Kq = 1 is our 'grain 
unit' . 

The results obtained from our model (Hig. [2j do agree with 
the analysis done on eight main financial indices, as shown 
in Fig. [3] Here, the time-series of the logarithmic returns 
(Fig.|3^) must first be mapped in an series of events as shown 
in Fig. [3j?. The non-Gaussian distribution of the logarithmic 
returns (Fig. [3};) were extracted from the logarithm returns of 
the original series of each index, similarly to what is done in 
Ref. lUfl. The characteristic heavy tails observed in |Q] are ob- 
served for short time lag (hours or smaller), wheres in Fig. [3}; 
the daily closure values are considered. The power-law be- 
havior of the avalanche size (Fig. [3}l) is indeed similar to the 
simulated results. Moreover the exponents m have all approx- 
imate values, plotted in Fig. [3^, around the simulated value 
m = 2.51 (solid line). 

This exponent can be derived analytically as follows. Being 
r the number of nodes in the multiplicative process that repre- 
sents the avalanche and taking the probability Pb r for a node 
to be part of it one can write r = LKo(ujki n )) 1 where L is 
the total number of nodes in the system. From r one gets the 
total number k jri of connections involved in the avalanche: 
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From Otter's theorem 11 111 for branching processes the dis- 
tribution for r is P(r) oc r~ 2 which leads to a degree distri- 
bution P(k) oc fc~2T and therefore the avalanche size distri- 
bution reads 



P(k > s) oc 



+00 



k 21 ' dk oc s 2 
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Equation dH) relates the exponent characterizing the network 
topology with the exponent taken from the avalanches, us- 
ing first principles in economy theory. Since, typically 7 6 
[2.1, 2.7] one arrives to the values of m observed in Fig. (f3]) 
and in our model (Fig. EJ;). 

Figure [4] shows how the overall index dynamic emerges 
from the underlying network mechanics. As agents connect 
each other by preferential attachment, the topology of the sys- 
tem is pushed to a power law degree distribution. On the other 
hand, avalanches push the system away from it. Thus, the sys- 
tem undergoes a structural fluctuation that generates a fat tail 
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indices can be explained by the different social-economic re- 
alities beneath them, dispite the fact that, in general, they all 
have the same behaviour. Our numerical model was based on 
the same principles to simulate drops on energy and repro- 
duces the drops on the financial indices with an error smaller 
than the differences between real indices. Interestingly, the 
exponent m = 2.5 appears in other contexts and physical 
systems, e.g. in Fiber bundle models II 1911 . which raises the 
question if the financial systems belong to some broader uni- 
versality class. Further, while the local crashes appear ran- 
domly in time, one can argue that with such a model predic- 
tions of broad avalanches and economic crisis may be possi- 
ble, similarly to what was recently demonstrated in granular 
systems 1 1011 . 



FIG. 4: (a) Time series of the total internal energy Ut of the agents. 
For the four instantes t\, t-x, ta and ti we show the observed cummu- 
lative degree distribution P(k* > k) oc k~ 1+1 yielding exponents 
(b) 7 = 2.14, (c) 7 = 2.25, (d) 7 = 2.30, (e) 7 = 2.11, all of 
them according to the theoretical prediction (check Eq. l[8}). The de- 
viations from the power-law for large k are due to the avalanches 
(crisis) in the system (see text). 

distribution of the index, expressed by Eq. ([8]l and shown in 
Fig. |3] rather then a gaussian one. Figure^ shows a typical 
set of successive Ut values taken from our model. The cum- 
mulative degree distribution P(k) at the marked instants t\-t^ 
are shown in Fig.[4j3-|4h. The dashed lines guide the eye for the 
scaling behavior observed at low degree. For large degree k 
the distribution deviates from the power-law, due to the drops 
of connections for agents experiencing an economical crash. 
Nonetheless, in every case the slope —7 + 1 of the dashed line 
yields values in the predicted range, namely [2.1,2.7]. 

We have showed that, based only on first principles of eco- 
nomic theory and assuming that agents form an open system 
of economic connections organized by preferential attachment 
mechanisms, one reaches the distribution of drops observed 
in financial markets indices, including stocks and interest rate 
options. Assuming that the preferential attachment mecha- 
nism is part of the growing of economic connections, the re- 
sulting self-similar topology allow us to assume that all econ- 
omy has the same topological structure as the sub-economy 
around financial markets and, thus, market indices are good 
proxies for the all economy. Since the exponent of the degree 
distribution is reflected directly on the exponent of the drop 
magnitude distribution, the diferences between the different 
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